The Astrophysical Journal, 682: 373, 2008 July 20 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 
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ABSTRACT 

Energetic flares are observed in the Galactic supermassive black hole Sagittarius A* from radio 
to X-ray wavelengths. On a few occasions, simultaneous flares have been detected in IR and X-ray 
observations, but clear counterparts at longer wavelengths have not been seen. We present a flare ob- 
served over several hours on 2006 July 17 with the Chandra X-Ray Observatory, the Keck II telescope, 
the Caltech Submillimeter Observatory, and the Submillimeter Array. All telescopes observed strong 
flare events, but the submillimeter peak is found to occur nearly 100 minutes after the X-ray peak. 
Submillimeter polarization data show linear polarization in the excess flare emission, increasing from 
9% to 17% as the flare passes through its peak, consistent with a transition from optically thick to 
thin synchrotron emission. The temporal and spectral behavior of the flare require that the energetic 
electrons responsible for the emission cool faster than expected from their radiative output. This 
is consistent with adiabatic cooling in an expanding emission region, with X-rays produced through 
self-Compton scattering, although not consistent with the simplest model of such expansion. We also 
present a submillimeter flare that followed a bright IR flare on 2005 July 31. Compared to 2006, 
this event had a larger peak IR flux and similar submillimeter flux, but it lacked measurable X-ray 
emission. It also showed a shorter delay between the IR and submillimeter peaks. Based on these 
events we propose a synchrotron and self-Compton model to relate the submillimeter lag and the 
variable IR/X-ray luminosity ratio. 

Subject headings: Galaxy: center - black hole physics - polarization 
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1. INTRODUCTION 

The radio, IR, and X-ray source Sagittarius A* is as- 
sociated with a supermassive black h ole at the center 
of our Galaxy (|Melia fc Falckd l2001h . Spectral mea- 
surements at all wavelengths where Sgr A* is not hid- 
den by confusion or Galactic absorption show it to be 
extremely underluminous for its mass, radiating just 
10~^LEdd- A variety of physical models have been 
shown to ade quately reproduce the quiescent spec trum 
of Sgr A* (e.g. lFalcke fc Markofl[2000t iMeha et al.ll2001l : 
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lYuan et al.ll2003 ^. Discrimination between the proposed 
accretion and outflow models will require information 
complementary to the spectral data. 

Since the discovery o f X-ray and IR flar e s in Sgr A* 
llBaganoff et al.l l2001l : iGenzel et all l2003t iGhez et al.l 
l2004f ). transient events have been studied extensively. 
Such observations have found that Sgr A* is highly vari- 
able, with increases in X-ray lumi nosity of up to 160 
times over the quiescent emission (iPorauet et al.ll2003D 
and smaller flares on hour timescales at longer wave- 
lengths (e.g. IGhez et al.l 12004 iMauerhan et al.l 120051 : 
iHerrnstein et al.ll2004f ). Because of the rapid modula- 
tion observed in the flaring emission, these events likely 
occur just outside the event horizon and may provide 
insight into the structure and conditions in the inner ac- 
cretion regions. Models for the flares have considered 
various mechanisms for injecting energy into the elec- 
trons, including stochast ic acceleration, shocks, and mag- 
netic reconnection (e.g. iMarkoff et al]|2001l : lYuan et al.l 
12004 iLiu et al.l l2006b[ ). The radiative processes re- 
sponsible for the flares at each wavelength have also 
been debated, leaving synchrotron and synchrotron self- 
Compton (SSC) emission as the most likely candidates 
for the IR and X-ray emission. 

Constraints on the emission processes have im- 
proved as more flares have been observed in the IR 
and X-ray bands. However, a great deal of un- 
certainty was generated by the conflicting measure- 
ments of the IR spectral index during the flares 
(IGhez et al.ll2005l: lEisenhauer et al]|2005l : lGillessen et all 
120061 : iKrabbe et al.ll2006 h. In particular, some previous 
attempts to explain the IR and X-ray spectra with syn- 
chrotron and SSC, respectively, have been forced to in- 
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elude complieations in order to explai n correlated vari- 
ations of the flux and spectrum (e.g. iLiu et alJ l2006"at 
lYusef-Zadeh et al.ll2006bt iBittner et al.l 120071 ). In a re- 
cent paper, Hornstein et al. (2007; hereafter H07) have 
used multi-band IR observations of several flares to show 
that, after corrections for stellar contamination, the spec- 
tral index of Sgr A* is roughly constant within and be- 
tween flares, with Si, oc v~'^'^ . This can be understood as 
optically thin synchrotron emission from a population of 
power-law electrons with an N{E) cx energy spec- 

trum. As discussed below, their findings can be used to 
make a strong case for SSC production of X-ray flares, 
as has been suggested by many authors. 

Previous considerations of flare emission have largely 
avoided the temporal evolution of the flares. An 
exception is the e xpandi ng plasmon model used by 
lYusef-Zadeh et ahl ()2006aD to explain delays between 
flares observed in two centimeter wavelength bands; the 
extension of this model to shorter wavelengths is dis- 
cussed in this paper. While understanding the flare cre- 
ation mechanism (without regard to the flare evolution) 
is an important goal on its own, much of the poten- 
tial of the flare measurements to constrain the struc- 
ture of Sgr A* comes from modeling the flare evolution 
in time and wavelength as the energized electrons cool 
and expand through the source. Initially, X-ray and IR 
flares lacked complementary information at other wave- 
lengths, limiting time- and frequency-domain studies to 
the information encoded in a single narrow band. To 
date, a handful of events have been detected simultanc- 
ously in X-rays and the IR (jEckart et al. 2004, 2006b; 
lYusef-Zadeh et aLll2006b[ ). yet because the IR and X-ray 
flares are observed to be simultaneous the flare evolution 
has received less attention than the peak spectra. 

A few flares have provided evidence for decaying 
millimeter and submillimctcr emiss ion following short- 
wavelength flares (IZhao et al.l l2004l: lEckart et al] l2006bl: 
lYusef-Zadeh et al.l l2006bir ~ Coordinated observations 
from centimeter to X-ray wavelengths are now routinely 
attempted to search for clear flare counterparts across as 
broad a wavelength range as possible. Here we present 
the first observations of a flare of Sgr A* detected at sub- 
millimeter, IR, and X-ray wavelengths. Using an array 
of telescopes (^ we are able to measure the amplitude, 
spectral index, and temporal structure of the flare in each 
band (Sj3]). We also report a second IR/submillimeter 
flare, detected in the same monitoring campaign, that 
lacks an X-ray counterpart. We find large delays between 
the time of the IR and X-ray flares and the submillimeter 
flares. In 21 attempt to constrain the emission pro- 
cesses and dynamics responsible for these and other flares 
observed in Sgr A*. We find that the timing, spectra, 
and energetics of the flares imply a synchrotron origin 
for the IR emission and a SSC X-ray generation mech- 
anism. The decay of these and the submillimeter flares 
also suggests that non-radiative cooling processes, such 
as adiabatic expansion, are essential. As an initial step 
toward understanding the structure of Sgr A* through 
the flare changes, we compare the present flare to an 
existing expansion model. Finally, we use simple scaling 
arguments to predict the relationship between the IR and 
X-ray flare luminosities and the delay of the submillime- 
ter counterpart. 

Throughout this paper we refer to spectral indices 



(a) using t he convention 5*^ cx v'^ . We assume the 
iReidI ()1993D distance to Sgr A*, 8 k pc, which is con- 
sistent with more rece nt results (e.g. iGhez et al.l [20031 : 
lEisenhauer et al.ll2003f) . 

2. OBSERVATIONS AND REDUCTION 

The data presented here were obtained as part of a 
2005-2006 campaign to monitor Sgr A* simultaneously 
across a broad range of wavelengths - these results en- 
compass data from four observatories spanning seven 
decades in wavelength. We report on two strong flares 
observed at IR and submillimeter wavelengths, only one 
of which was accompanied by an X-ray flare. The tempo- 
ral coverage at the various observatories is shown in Fig- 
ure [TJ Details of the individual observations and analysis 
techniques are discussed in the following sections. 

2.1. X-ray Data 

Th e Chandra X-Ray Observatory ([Weisskopf et al.l 
Il996l ) observed the Galactic center on both 2005 July 
30/31 and 2006 July 17 using the ACIS imaging array 
(jCarmire et al.|[2003f ). The observations were timed to 
span the window of Sgr A* visibility from Mauna Kea for 
coordination with telescopes there. Obser vational details 
and analysis procedures followed those of lBaganoff et al.l 
(|2001ll2003D . In particular, Sgr A* photometry was ob- 
tained from 2—8 keV counts within 1.5", after subtrac- 
tion of a background derived from a 2—4" annulus with 
point sources and structures excluded. 

2.2. IR Data 

The W. M. Keck II 10 m telescope observed the Galac- 
tic center using the NIRC2 (PI: K. Matthews) near-IR 
camera and the la.ser gu i de star adaptive opti cs system 
(jWizinowich et all 120001 : Ivan Dam etaH l2006f ) on 2006 
July 17. Observations were alternately made in the K' 
(Ao = 2.12, AA = 0.35 /im) and L' (Ao = 3.78, AA=0.70 
fim) photometric bands, with exposure times of 28 and 
30 s, respectively, each cycle. The observations span 
187 minutes, with 16 minutes of data lost to instrument 
problems. Additional data were obtained on 2005 July 31 
in the H (Ao = 1.63, AA=0.30 /im), K' , and L' bands, cy- 
cling through 22.2, 28, and 30 s exposures in these bands 
every 3 minutes. Within the 113 minutes of observations, 
9 minutes were lost to telescope problems. The dead time 
between frames on these nights was typically less than 1 
minute. Seeing on both nights was excellent; the reso- 
lution achieved at H and K' was 62—65 mas (FWHM) 
and 80—82 mas at L'. We refer the reader to iII07l for 
additional details. 

2.3. Submillimeter Data 

The Caltech Submillimeter Observatory (CSO) 
SHARC-II observation and a nalysi s methods are de- 
scribed bv lYusef- Zadeh et al.l ()2008[ ). with attention to 
the 2006 July 17 observations at 850 /im. Observations 
were also made on 2005 July 31 at 350, 450, and 
850^lm, for which the CSO has 8.5", 10", and 20" 
resolution, respectively. For the 2005 observations, 
850 fim calibration was derived from Callisto (10.3 Jy) 
and Neptune (27.7 Jy), with an estimated uncertainty of 
10%. Confusion causes an additional ~1 Jy uncertainty 
in the absolute flux density of Sgr A*. At 450 /im. 
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Arp 220 (6.3 Jy), Callisto (35 Jy) and Neptune (67 Jy) 
were used for absolute calibration, with an estimated 
uncertainty of 25%. Confusion causes an additional 
~0.5 Jy uncertainty in the absolute flux density of 
Sgr A* at 450 /xm. At 350 ^iul, Arp 220 (10 Jy) and 
Neptune (93 Jy) were used for absolute calibration, 
with an estimated uncertainty of 25% and a confusion 
uncertainty of ~1 Jy. 

Submillimeter Array (SMA) observations of Sgr A* 
were made on 2005 July 31 and 2006 July 17 (UT). 
In 2005, seven antennas were used in the SMA "com- 
pact north" configuration at 1.32 mm wavelength 
(226.9 GHz), covering baseline lengths of 5 — 53kA and 
yielding a synthesized beam of 3.8"x2.1"after a 6 hr 
track. During the track the zenith opacity varied between 
0.05 and 0.08. The SMA polarimetry system (jMarrond 
[20061) was installed for these observations in order to con- 
vert the linearly polarized SMA feeds to circular polar- 
ization sensitivity. This removes the possibility of con- 
fusing linear polarization modulation with total intensity 
variations. Gain calibration was derived from the quasar 
J1733-130, while J1744-312, just 2.3° from Sgr A*, was 
used as a comparison source to verify the calibration. 
The flux density scale was determined from Uranus, with 
an uncertainty of 15%. In 2006, seven antennas were used 
in the "very-extended" configuration, yielding baselines 
of 27— 390kA and a synthesized beam of 0.6"x0.5"after 
a 6.5 hr track on Sgr A*. The observing wavelength 
was the same as in 2005, while the zenith opacity was 
0.10. For these data the polarimetry system was used 
to make full polarizatio n measurements according to the 
procedures described in M arrone et al.l (|200 6a!) . In order 
to sample all cross-correlations of left and right circular 
polarization on all baselines the feed polarizations were 
modulated in a coordinated pattern with a 4 minute du- 
ration; this cycle time set the minimum length of the po- 
larization samples. Instrumental polarization calibration 
was obtained through observations of the quasar 3C 279, 
yielding measurements consistent with those obtained in 
previ ous observ ations at this frequency. As described 
in M arrone et a l. (2007), the calibration precision limits 
false linear polarization signals to 0.2%. Gain calibra- 
tion was derived from J1626-298 and J1924-292, with 
J1733— 130 as a verification source. Callisto was used for 
absolute calibration, with an uncertainty of 15%. For 
both epochs, the complex calibrator gains were applied 
to the Sgr A* data and then Sgr A* was used for phase- 
only self-calibration. Projected baselines shorter than 
20 kA were excluded from this procedure because of con- 
tamination from extended emission around Sgr A*. Flux 
density measurements were obtained for each time inter- 
val (4 minute on-source cycles) by fitting a point source 
to the calibrated visibilities. Flux density uncertainties 
were adjusted to account for the precision of the cal- 
ibrator gain measurements, while the overall flux den- 
sity scale uncertainties reported above were not included 
because they should be common to all time intervals. 
Figure [2] shows the SMA light curves for both epochs, 
including calibration and verification sources. 

3. RESULTS 

3.1. Flare Amplitude and Duration 

Figure [3] shows the light curve observed at submillime- 
ter, IR, and X-ray wavelengths on 2006 July 17. All three 



bands (four telescopes) show a fiare between 6 and 8 hr 
UT. Assuming that the events seen at these wavelengths 
are related, this is the first flare of Sgr A* to be observed 
in all of these bands. 

The X-ray flare, centered around 06:10 UT, has a 
FWHM of 31 minutes and a FWZP of roughly 1 hr. 
At its peak, this flare has a 2 — 8 keV luminosity of 
4.0 X lO^** ergs s~^, approximately 20 times the quiescent 
X-ray luminosity of Sgr A*. The integrated emission of 
the flare has a spectral index of a = O.O^J g [photon index 
of r = 1.0 for N{E) = E~^], implying a monochromatic 
luminosity {vL^) of 7 L© at 4 keV. Flares of this ampli- 
tude or larger have been observed on six occasions in the 
past jB aganoffet al. 2001, 2002; Goldwurm et al. 2003|; 
iPorquet et al.ll2003HBelanger et al.ll2005l ). corresponding 
to a rate of around 0.6 ± 0.3 day~^. 

The IR observations begin 36 minutes after the peak 
of the X-ray fiare and Sgr A* is initially a factor of a 
few brighter than the minimum emission observed over 
the night. The 7 (7.5) mJy peak observed at K' {L') 
corresponds to 20 (12) Lq {vL^). Within 50 minutes, 85 
minutes after the X-ray peak, the emission decays to a 
low level (2-3 mJy). Throughout the IR flare the K'-L' 
spectral index is approximately —0.51 ib 0.14, consistent 
with other IR flares discussed in |H07| . The spectral in- 
dex between the K' and X-ray peaks is —1.21, although 
more negative indices are allowed because the IR peak 
may have been significantly brighter. Assuming that the 
K' peak was comparable to the largest flares observed 
to date, ~ 12 mJy, the spectral index would be —1.28. 
In previo us observatio n s of X- ray fl ares with IR coun- 
terpart s, lEckart et al.l (|2006b( ) and lYusef-Zadeh et all 
(I2006bf) found ax-x = -1.12 and an-x = —1 .3, re - 
spectively. Assuming an IR spectral index of —0.6 ()H07[ ). 
the latter is equivalent to aK'-x = —1.2. 

Although the submillimeter observations span the X- 
ray and IR flares, there is no submillimeter flare apparent 
at the time of the maxima in these bands. Prior to the 
X-ray fiare, both telescopes show a small (0.2 Jy, < 10% 
fractional change) rise and fall in fiux density. Due to an 
unfortunate coincidence neither telescope was observing 
Sgr A* precisely at the peak of the X-ray flare, but there 
is no suggestion of a missed increase in emission from 
the data immediately before or after the gap. However, 
a large (1 Jy) fiare is seen at both wavelengths, peaking 
more than an hour after the X-ray fiare. At 1.3 mm and 
850 /xm the monochromatic luminosities of 1 Jy flares are 
4.6 and 7.0 Lq, respectively. Events of this magnitude 
have been se en in previous observations at 1.3 mm and 
850 /xm (e.g. [M arrone et al.ll2006i lEckart et al.ll2006bl : 
lYusef-Zadeh eFa l. _2006J3); they occur with a frequency 
of 1.2 day~^ based on 20 epochs since 2004. The decay 
of this fiare is well-approximated at both wavelengths by 
an exponential with a time constant of 2 hr. A similar de- 
cay was also suggested by the 850 /im data presented by 
lEckart et al.l (12P06b). The spectral index of the fiaring 
component is tough to determine because of the absolute 
calibration uncertainty and the difficulty in determining 
the non-flaring fiux. Assuming that the minimum flux 
density observed at each wavelength represents the stable 
component, the submillimeter spectral index during the 
flare rise (07:00-07:30 UT) is asubmm = -0.1 ±0.2 ±0.4, 
with a mean of 0.4 ± 0.1 ± 0.4 after the flare peak. For 
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each spectral index we separate the errors resuhing from 
the measurement error (first) from the constant error due 
to uncertainty in the absolute calibration of the two ob- 
servations (second). The change in spectral index across 
the peak of the flare is an increase of 0.5 ± 0.2. 

On 2005 July 31 we also observed a strong IR flare, 
among the brightest yet detected (Figured]). It was ac- 
companied by a 1.3 mm flare of similar amplitude to 
that of 2006 July 17. The IR and submilhmeter-t o-IR 
indices are very similar to those in the 2006 flare; iHOTi 
report a spectral index of —0.62 ± 0.21 between K' and 
L'an d ai,'imm-K' — —0.7 in both epochs. As noted by 
IH07L there is no appreciable change in X-ray flux dur- 
ing these observations, despite coverage beginning more 
than 10 hr before the start of the Keck IR data. The 
non-detection of X-ray emission places an upper limit of 
uk'-x < —1.50. 

3.2. Correlation Analysis 

The peak of the 2006 July 17 X-ray flare occurred be- 
fore the beginning of our IR observations. The probabil- 
ity of the IR and X-ray flares coinciding by chance within 
this time interval is non-ne gligib le given the observed IR 
flare rate, as discussed bv IH07| . However, all previous 
X-ray flares that have occurred during IR observations 
have been accompanied by an IR flare, with no mea- 
surable time del ay between the two wavelength bands 
(< 10 minutes: lEckart et all [2006bl: lYusef-Zadeh et all 
I2006bt lEckart et all 120081 ). The apparent flare peak 
at the beginning of the IR observations is consistent 
with substructure o bserved in previous IR events (e.g. 
lEckart et al.]|2006bO . We therefore assume that the IR 
and X-ray peaks are coincident and expect that the max- 
imum IR flux density was greater than the ~7 mJy at 
the start of these observations. We refer to lH07l for fur- 
ther discussion. 

Neither of the flare events in 2005 and 2006, as marked 
by the IR and X-ray emission, shows coincident sub- 
millimeter activity. Both, however, show submillime- 
ter flares of unusual amplitude after the X-ray or IR 
emission peak. The apparent delay between the sub- 
millimeter and IR/X-ray flares makes the assertion of a 
relationship between these events even more uncertain 
than the IR/X-ray connection described above, but cir- 
cumstantial evidence of a relationship is building. From 
campaigns between 2004 and 2006 there are approxi- 
mately 52 hr of joint X-ray/submillimeter observations 
of Sgr A* yieldi ng just one X-ray f l are, the 2006 flare 
presen ted here (jEckart et al.l l2006bl : lYusef-Zadeh et all 
I2006bl: F. K. Baganoff et al., in preparation). A 2004 
flare that occurred 2.3 hr before the start of the sub- 
millimeter observations was also followed by a 0.8 Jy 
decline in 870 /im flux over the firs t 2 hr of the sub- 
millimeter light curve (jEckart et al.l l2006bD . A simi- 
lar number of hours of simultaneou s IR /submill imeter 
measurements (jEckart et al.l l2006bl : l^scf-Zadch et al.l 
I2006U this work) have produced three instances of IR 
flares followed by submillimeter flares. In the case of 
the 2005 flare presented her e and t he 2004 September 
4 flare in lYusef-Zadeh et ahl (j2006b[ ) , the submillimeter 
event occurs after a large IR flare but precedes a smaller 
flare. In these cases it is not clear which IR event to as- 
sociate with the submillimeter, if any, but we note that 
in all five of these cases the submillimeter flare follows 



the largest event observed at the shorter wavelengths. If 
the X-ray/IR events are unrelated to the submillimeter 
we would expect an equal number of flares before and 
after the X-ray/IR flares. We therefore proceed on the 
assumption that the two submillimeter flares presented 
here are related to the X-ray/IR flares. 

In Figure |5| we show the cross-correlation of the 2006 
submillimeter and X-ray light curves. We have employed 
the z-tran sform discr e te co rrelation function (ZDCF) 
analysis of lAlexandeil (jl997f ) in order to treat properly 
the irregular sampling of these data sets. We find no 
significant delay between the 1.3 mm and 850 /xm light 
curves, with the 850 /zm peak leading by 2 ± 12 minutes. 
Cross-correlation with the X-ray light curve indicates de- 
lays of 96 ± 14 and 97 ± 17 minutes for the 850 /im and 
1.3 mm data. The cross-correlation of the 2005 IR and 
1.3 mm data is also shown in Figure (top), where the 
IR flux is represented by the spectral-average light curve, 
obtained by scaling the H and K' flux densities to the L' 
band through the factor {v/vli)'^, whe re a = —0.62, the 
mean K'-V spectral index of the flare (jHOTT ) . This com- 
posite light curve leads the 1.3 mm flare by 20 ± 5 min- 
utes; cross correlation with each individual IR light curve 
yields similar results and the inter-correlations of the IR 
light curves show no evidence for relative delays. The 
second peak in the cross correlation is spurious, arising 
from the chance alignment of the gap in the SMA data 
with a minimum in the IR light curve. 

The lag between the submillimeter and X-ray flares 
in the 2006 event is nearly 80 minutes longer than the 
submillimctcr-IR lag in the 2005 flare. However, because 
the 2005 flare shows no X-ray emission and we lack IR 
coverage at the expected peak of the 2006 flare, we can- 
not compare cross-correlations of the same pair of wave- 
lengths between the two flares. If the plateau at the 
beginning of the 2006 IR data truly represents the peak 
of the IR flare, the delay between IR and submillimeter 
would be 40—45 minutes shorter, although still measur- 
ably longer than that observed in 2005. 

3.3. Flare Polarization 

Although linear polarization has been detected in 
Sgr A* at submilli m eter and near-IR wavelengths (e.g. 
Aitken et al] 120001: iBower et al.l l2003t iMarrone et all 
2006at lEckart et al.ll2006al ). of the observations presented 
here only the 2006 July 17 SMA observations were de- 
signed to measure polarization. The 1.3 mm polarization 
light curve for that epoch is shown in Figure |3| (6otto77i). 
The fractional polarization varies from 1%— 2% at the 
start of the track to as much as 8%— 9%. The polariza- 
tion position angle varies between 90° and 130° in the 
4 minute samples. 

If the submillimeter flare emission arises from the syn- 
chrotron process, the flare might be expected to be highly 
polarized. To examine the flare polarization, we re-bin 
the data in half-hour intervals (typically, four 4 minute 
observing cycles) and subtract the total intensity (/o) 
and polarization (Qoj Uq) averaged over the four sam- 
ples that precede the onset of the submillimeter flare. 
Although Sgr A* often shows dramatic polarization mod- 
ulation (magni tude and dire ction) even during periods 
of quiescence (jMarrond 12006 ) , making the assumption 
of a single Qo and possibly unreliable, the result- 
ing background subtracted light curve (Figure [6]) reveals 
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interesting changes. As the excess Stokes / rises and 
falls a polarization component (Pexcess) also appears and 
fades, suggesting that the flare emission is significantly 
polarized. Previous IR and centimeter-wave observa- 
tions of Sgr A * have also shown evidence of polarized 
flare emission (lEckart et al.ll2006at iMever et al1 l2006bl: 
iTrippe et al] 120071: 1 Yusef-Zadeh et al.ll2007f) . The polar- 
ization fraction of the excess emission (niexcess) is ob- 
served to increase from 9.4 ± 1.9% while the flare inten- 
sity is increasing (the first two bins after the flare onset) 
to a weighted average of 16.5 ± 2.3% after the peak (ex- 
clusion of the last bin causes an insignificant change in 
this average). This increase is consistent with a syn- 
chrotron flare that is evolving from optically thick to 
optically thin, assuming a power-law electron distribu- 
tion with N{E) cx E^P and p > -0.45. For the elec- 
tron index indicated by the constant IR spectral index, 
p = 2.2, the polarization fraction would be expected to 
change from 11% to 71% through this transition if the 
flaring region lacked any appreciable random magnetic 
field component. The smaller change observed here sug- 
gests that there is significant disorder in the field revealed 
as the flare becomes optically thin, or substantial internal 
Faraday rotation at 1.3 mm. 

The variation of the excess Q and U through the flare 
represents a rotation of the polarization, the total ex- 
cess polarization (Figurc[6l bottom panel, circles) remains 
nearly constant. Comparing the data point on the rising 
edge of the flare with the six after the peak we find that 
the polarization angle changes by 40°, not as large as the 
expected 90° change through a transition from optically 
thick to thin synchrotron emission. However, the magni- 
tude of this change depends strongly on the choice of Qo 
and Uo and could be made to agree with the prediction if 
these quantities are slightly more negative than assumed. 

4. DISCUSSION 

4.1. Emission Mechanisms and Electron Cooling 

After several years of coordinated multi-wavelength 
monitoring of Sgr A*, the physical conditions and mech- 
anisms responsible for its flaring in are becoming clear. 
The IR observations of Sgr A* in flaring and quiescent 
states by iH07i show a consistent am — —0.6 spectrum, 
independent of the instantaneous flux density and its 
derivative. The spectral index suggests that the IR pho- 
tons are optically thin synchrotron emission from power- 
law electrons [N{E) cx E'P] with p ^ (1 - 2a) = 2.2. 
Moreover, the stability of the spectral index as the flares 
decay is inconsistent with the Aa > —1/2 change ex- 
pected if the decay results fro m radiative cooling of 
the electrons (jPacholcz vkl 1 1 9 70() . The ele ctron cooling 
timescale due to synchrotron losses is (e.g. iKrolikI 0*9991 ) 



tcool = 1.3 X lO^^iy-^/^ 



B-^^\ (1) 

where the frequency (i') is in Hz and the magnetic field 
(B) is in G. Assuming that after the flare peak the IR- 
emitting electrons are no longer produced in large num- 
bers and can no longer hide a change in spectrum, the 
25 minute decay of the IR flares limits the field in the 
emission region to ~20 G. At this field strength, electrons 
emitting at K' have a Lorentz factor of 7 > 1600. 

Measurements of bright radio and submillimeter flares 
also imply that we are observing synchrotron flares that 



decay due to nonradiative electron cooling. The strongly 
polarized flare emission shown in Figure[3]is suggestive of 
a synchrotron origin. In both of the submillimeter flares 
considered here the excess flux fades within 2 hr, much 
more quickly than could be explained by synchrotron 
losses. Equation ([T]) predicts that submillimeter-emitting 
electrons should cool 20 times more slowly than those ob- 
served in the IR bands, very different from the observed 
factor of a few difference in decay time. Similarly rapid 
decay has been observed in flares at lower frequencies 
()Yusef-Zadeh et al.l [2006al ). The behavior of the long- 
wavelength flares and the achromaticity of the IR de- 
cay imply that an energy-independent process, such as 
expansion, dominates the energy loss. Magnetic flux- 
conserving expansion also reduces the magnetic field and 
therefore could allow a somewhat higher initial B and 
smaller Lorentz factor than those quoted above. 

Between 2000 and 2006, Chandra and XMM-Newton 
have found 11 significa nt increases i n the X-ray 
luminosity of Sgr A* dB aganoff ct al.' '200T, '2OOI 
Goldwurm et al.ll2003t jPorcmet et al. 2003 : ^Eckart et all 
2OO5I: lEckart et al] l2006bt this 
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work). Typically these flares last for 0.5—2 hr, much 
longer than the synchrotron lifetime for reasonable es- 
timates of the magnetic field strength, so production 
of the X-ray flares through direct synchrotron emission 
would require sustained injection of high-energy elec- 
trons throughout t he flare (e.g. iBaganoff et al.l l2001l : 
iMarkoff et all I200I . On every occasion where IR 
data have been available, IR coun terparts to these 
flares have been observ ed (|Eckart e t al. 2004, 20061^; 
lYusef-Zadeh et al]|2006bl : this work). For the two flares 
with the best data, those where the flare rise and fall 
was observed in both bands, there is no significant de- 
lay between the two wavelengths (Eckar t et al.l l2006bl : 
[Yuscf-Zadeh ct al. 2006b). H07 also noted the correspon- 
dence between the X-ray spectral indices {ax ) and their 
mean am for all but the brightest X-ray flare observed 
to date. Finally, as discussed in ii3.1l the spectral index 
between IR and X-ray wavelengths is variable [uk-x 
ranges from —1.1 to < —1.5) but is reliably more neg- 
ative than the spectral indices within the IR or X-ray 
bands. When taken together, these points demonstrate 
that the X-rays are produced through inverse-Compton 
scattering of the lower frequency spectrum (see also H07; 
IXusef-Zadeh et al. 2006b). 

4.2. Expanding Plasmon Flare Evolution Model 

lYusef-Zadeh et all (|2006al ) proposed that the tem- 
poral and spectral behavior of centimeter-wavelength 
flares in Sgr A* could be explained in an ex pand- 
ing synchrotron plasmon pictu re (Shklov sklil Il960l : 
iPaulinv-Toth fc Kellerniannlll966f). following the fo rmu- 
lation of Ivan der LaanI ( 19661 ). lEckart et al.l ()2006b[ ) also 
proposed an expansion model, although with a less di- 
rect connection to previous work. Fundamental to this 
model is the adiabatic cooling of electrons in the plasmon 
and the flux-conserving diminution of the magnetic field, 
which provide the nonradiative decreases in synchrotron 
output that we require. The model predicts smaller and 
later flare peaks at longer wavelengths, with the spectral 
indices characteristic of, respectively, optically thick and 
thin synchrotron [a = 2.5 and (1 — p)/2] before and after 
the flare peak at a given wavelength. 
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This particular model can be tested in new ways with 
the 2006 flare because we have observed the flare at 
two optically thick wavelengths (1.3 mm and 850 ^m), 
know the electron spectral index from the IR observa- 
tions (p — 2.2), and from the X-ray data can pinpoint the 
time at which the putative expansion was initiated. We 
found above that the submillimeter spectral index prior 
to the flare peak is —0.1 ± 0.5 in the flaring component, 
inconsistent with the expected value of 2.5. This latter 
number is a direct result of the assumption of a homo- 
geneous plasmon, but allowing variations in the electron 
density and magnetic field with optical depth, as in a je t 
or other inhomogeneous structure (e.g. Ide BruvnlfT"976l ). 
is well known to produce arbitrary spectral shapes. Sim- 
ilarly, the optically thin spectral index was found to be 
0.4 ± 0.5, just marginally consiste nt with the —0.6 ex- 
pected from the electron spectrum. iDentI (|1968l ) pointed 
out that the light curve maxima at two wavelengths sat- 

isfy Sjn,i/Sm,2 = {iyi/^2)\ where S^{7p + 3)/ (4p + 6) 
(ranging from 1 to 1.46 for p = 1—5). For our submil- 
limeter data and p = 2.2, we expect the 850 ^m peak to 
be 1.7 times brighter than the 1.3 mm peak. We instead 
find the amplitude of the flares in these two bands to 
be very similar, (S'gsoMm/'S'i.Smm) = 1.15 ± 0.15, consis- 
tent with p ~ 0. Finally, within this expansion model 
the relative timing of the flare peaks at these two wave- 
lengths is ih/h) - (1^1^)% e - -(p + 4)//3(4p + 6), 
for expansion as r ^ t^. Here the flare peaks at short 
wavelengths (IR/X-ray) at t = to, the scale time, mea- 
surable at soine opti cally thick wavelength as to = 3S/S 
()van der LaanI [T966D . Derivation of the scale time is 
quite uncertain due to the quiescent emission and short 
rise time, but from the 850 /im light curve we infer 
ti.3 mm— ^850 fim — 34 minutcs and setting to — places a 
lower limit of 18 minutes on the expected delay for linear 
expansion (/3 = 1). The latter is marginally inconsistent 
with the observed delay, while the former is discrepant at 
3 a. We also note that in this model the minimum delay 
between 7 mm and 850 /im should be 135 minutes, so in 
this context we do not expect any relation betwee n flares 
observed at 7 and 14 mm in lYusef-Zadeh"eraII (l2008h 
and the large 1.3 mm/850 /im flare considered here. 

Perhaps a more important problem is revealed by con- 
sidering the expansion r ate expected for t h e relat ivis- 
tic plasmon. Although lYusef-Zadeh et all (|2006aD in- 
voke an expansion speed of 0.02c, the sound speed near 
the black hole should approach the relativistic limit of 
c/V3. In the case of the 2006 flare, the submillime- 
ter peaks occur nearly 100 minutes after the event that 
initiated the putative expansion, implying an expansion 
distance of 10^'' cm (100 rs). Although the expansion 
speed may decrease as the plasmon entrains material, 
this estimate is nearly 2 orders of magnitude larger than 
the likely size of Sgr A* at submillimeter wavelengths. 
Extrap olations of rnillimeter-wavelength VLBI measure- 
ment s (jBower et al.ll2004l: IShen et al.ll2005t IBc 



jower 



et al.l 



l2006l ) suggest an intrinsic quiescent source size of ^^2 rs 
at 850 ^m. Further evidence for the small submillime- 
ter size comes from SED measurements that find the 
turnover in the submillimeter spectrum expected from 
the transi tion to opti cally thin emission (jMarrone et al.l 
[2006b; M arronel[2006l ). Therefore, 25% of the luminosity 
of Sgr A* near 1 mm would need to be produced by a 



plasmon that has roughly 50^ times the surface area of 
the quiescent source, implying a remarkably low bright- 
ness temperature in the plasmon. 

The expanded size can be transformed back to an ini- 
ti al size through the op acity law. Under the assumptions 
of lvan der LaanI (|1966f l. the synchrotron opacity depends 
on the frequency and expanded size as 

where tq is the opacity at a reference wavelength Aq and 
(R/Ro) is the expansion factor. Using p — 2.2, assum- 
ing that initially T3.8^,„ < 0.5 to match the spectral in- 
dex stability constraint ()H07l ) and that at the time of 
the 850 /im peak r = 1 .6 as predicted by the mo del for 
this electron spectrum (jYusef-Zadeh et ani2006af l. equa- 
tion ^ shows that (R/Rq) < 8.3. If the plasmon ex- 
pands by no more than this factor before reaching a size 
of 100 rg at the time of the 850 /im peak, the initial size 
is at least 12 rg. This source size would dramatically 
overproduce the observed IR luminosity unless the den- 
sity were very low (rig ~ 10^ cm^'^) or the assumption 
of homogeneity were removed. 

4.3. Other Dynamic Flare Models 

It is clear that although the Ivan der LaanI (|1966D plas- 
mon model grossly predicts some features observed in 
this flare, it is inadequate to describe the data presented 
here. However, this model is just one realization of a fam- 
ily of models that describe the scaling of energy and the 
magnetic field under expansion. The physics and geome- 
try/dimensionality of the expansio n may prescribe other 
scaling relatio ns fe.g. | Ko nigl"1981') or more complicated 
variations (e.g. lFalcke fc Markofl 2000). It is established 
above that the properties of Sgr A* flares require non- 
radiative electron cooling f !j4.ip . Detailed models that 
describe the density and field structure in the accretion 
flow or outflow can also predict the evolution of an ex- 
panding region, so time-resolved multi-band flare obser- 
vations can directly test the structure of these models. 

It has often been argued that the submillimeter spec- 
trum of Sgr A* is dominated by an electron component 
that is not signifi cant at other wavelen gths, the "sub- 
millimeter bump" (jMelia fc Falckdl2001l ). Our SMA and 
CSO observations fall on the long-wavelength side of the 
peak of this bump, where the synchrotron emission from 
this component is optically thick. It is therefore pos- 
sible for the properties of the submillimeter fiare to be 
significantly altered by the excess opacity. For exam- 
ple, for some period of time the ambient submillime- 
ter bump ele ctron population, often taken to be thermal 
([Yuan et al.ll2003l ). could enshroud otherwise observable 
emission from the flaring region. However, the submil- 
limeter photosphere is believed to be small (few rg) based 
on extrapolated VLBI size measurements and it there- 
fore seems unlikely that this mechanism can hide an ex- 
panding blob for long. It is also possible that the flare 
electrons that produce submillimeter radiation are not 
injected into a power-law tail but instead are heated into 
a thermal spectrum. In this case, the differing depen- 
dence of the thermal synchrotron absorption coefficient 
on the source properties will change the simple relation- 
ship between opacity and expansion derived for power- 
law electron distributions. We have made no attempt to 
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treat these possibilities, although they are likely to be 
very important for proper modeling of flares with sub- 
millimeter observations. 

Observations of repeated structures in IR and 
X-ray flares (e.g. lEckartet"an l2006at iMever et all 
I2006bt iBelanger et"all 120061 ) have often been at- 
tributed to plasma "hot spots" orbiting the black 
hole fe.g.lBroderick fc L oetj'200a iMarrone et an i2006bt 
IMever et al.l l2006at |Trippe ct alj I2007D . In these in- 
terpretations, intensity and polarization features with 
~20 minute cycle times are ascribed to orbital motion, 
with several cycles observed in some flares. If these fea- 
tures are to persist for multiple orbits they must not 
expand significantly. However, the decay timescales for 
X-ray flares, indicative of expansion, are typically com- 
parable to a single orbital period and conflict with the 
required plasmon confinement. Unless separate mech- 
anisms are invoked for the "periodic" single-band flares 
and multi- wavelength flares shown here and elsewhere, it 
is unlikely that such hot spots survive for several orbits. 

4.4. X-ray Emission and the Submillimeter Delay 

There are two striking differences between the 2005 
and 2006 flares in Figures [3] and 21 First, although the 
2005 IR flare reaches twice the peak (observed) flux of 
the 2006 flare, it shows no measurable X-ray emission 
above that from the quiescent extended component. Sec- 
ond, the delay between the short- wavelength and sub- 
millimeter flares in 2005 is much shorter than in 2006, 
although we cannot rule out that the submillimeter flare 
is related to the IR flare seen around 7UT rather than 
the much stronger flare at 8UT. Presuming that the X- 
rays arise from inverse-Compton processing of the longer 
wavelength spectrum and that the late appearance of 
the submillimeter emission results from optical depth 
changes, we use a simple synchrotron-SSC source model 
to estimate how the X-ray/IR ratio and submillimeter 
delay should be related. 

We begin with a homogeneous spherical synchrotron 
source of radius R, electron density n^, and magnetic 
fleld B. We assume a power-law distribution of electrons 
between 7mm and jmax, ^ (l) oc 7"^, with p = 2.2 
as determined from the IR spectrum. The scaling of 
the synchrotr on and SSC spectra of suc h a source were 
described by iBloom fc Marscheii (|1996f ). and we follow 
their analysis here. Approximating the spectrum of a 
single electron of energy 7 by a delta function at the 
characteristic emission frequency = 2.8j^B MHz, the 
synchrotron spectrum of the source at optically thin 
frequencies is then proportional to 

cx n,R^B<^'+Py^i^<^'~Py\ (3) 

The SSC spectrum {Sl'~^) is proportional to the Thomson 
optical depth of the sphere (roughly UeRdT) times , 

The ratio of these two equations, namely S^'-^ / oc rieR, 
provides an explanation for the variation in am-x noted 
in ^13. II differences in the density and size of the flaring 
region from flare to flare. This model preserves the spec- 
tral similarity of the X-ray and IR flares, matching the 
observations. 

We can numerically compare the observed IR and 
X-ray flare emission to our spherical source model by 



adapting the publi c ly ava ilable synchrotron-SSC code of 
iKrawczvnski et all ()2004f ). This code assumes a spherical 
emission region of radius R moving at an angle from 
the observer's line of sight at speed f3 = v/c, yielding the 
conventional Doppler parameter S = 1/ [7(1 — Pcos9)]. 
The electron spectrum is specified as a broken power law 
distribution and the density, magnetic field, and Doppler 
factor are also free. We have modified the code in small 
ways to suit our Galactic application, rather than the 
modeling of high-energy blazar spectra that led to its 
development. In Figure [7] we show three models that 
match the IR and X-ray spectra in the 2005 and 2006 
flares. The parameters for the models are shown in Ta- 
ble [TJ The model is underconstrained by the available 
data, so we hold the Doppler factor {6 = 1.8) and the 
range of electron energies (7^™ = 1, ^max = 3 x 10'') 
fixed. Models 1 and 2 match the 2006 fiare at its IR/X- 
ray peak (where the peak IR flux densities are assumed 
to match those of the 2005 flare) and near the start of 
the IR data, respectively. These differ by a 1.4 times 
adiabatic expansion, with the magnetic field strength 
held constant. Model 3 also reproduces the peak 2005 
IR flux densities but evades the X-ray upper limit be- 
cause of its smaller optical depth to Compton scattering. 
The model parameters are plausible for Sgr A*, although 
polarizati on measurements likely prefer smaller densi- 
ties (e.g. iMarrone et all l2007t iLoeb fc WaxmanI |2007^ . 
This toy model is driven to higher densities by the need 
to reproduce the X-ray emission entirely through self- 
Compton scattering; a more complete model of Sgr A* 
would include the quiescent submillimeter emission and 
these additional seed photons and scattering electrons 
would therefore permit a smaller flaring density. 

The delay between the flare times at optically thin 
and thick frequencies depends on the initial optical depth 
(■'"1^,0) at the frequency of interest and its rate of change. 
The synchrotron opacity o f the model sphere scales as 
(jRvbicki fc LightmanllTOTj) 

Calculation of a "lifetime" for this opacity requires that 
we introduce some relationship between the quantities 
in equation (O and time. The simplest procedure is to 
impose a pow er-law dependence of the radius on time, 
R oc t^, as in Ivan der LaanI (|1966[ ). with the magnetic 
field and density at fixed energy (or energy density) scal- 
ing as B oc R''" oc t^''^ and n oc i?''" oc t^^'^" . Inserting 
these equations into equation (O and its derivative, we 
find that the opacity decreases according to 



(6) 



In the case of the Ivan der LaanI (|1966f ) model, fc^ = — 2 
and fc„ = —2 — p, so /j, = — (3 + 2p). Then the time T^, 
required to reduce the source opacity to unity, the delay 
between the initial flare and the peak at frequency v, is 



(7) 

Combining equations ([5]) and ([7]), the delay depends on 
the initial parameters of the source according to 



T^oc UeflRoB^Q' 



(P+2)/2^-(4+j9)/2 
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i^^(P+2)/2,,-(4+p)/2 



We can use the observed X-ray and IR flux densities 
in place of Sl'^' and ; normalization factors that de- 
pend on frequency in equations ([3]) and ^ will cancel in 
comparisons between flares because the observing wave- 
lengths do not vary. 

The two flares presented here can be used to exam- 
ine the plausibility of this relationship. However, be- 
cause we can only compare these flares through ratios of 
their properties, we cannot test the model without ad- 
ditional observations of flares having submillimeter and 
X-ray and/or IR counterparts. Normalized to the qui- 
escent X-ray flux of Sgr A*, the 2006 flare represented 
a factor of 20 increase, while the 2 005 flare produced 
< 1.2 times the quiescent flux (|H07D . The observed IR 
{K') peak flux density was 12 mJy in 2005 and 7 mJy 
in 2006, but it is possible that the IR flux density was 
comparable to or even greater than 12 mJy in 2006 be- 
fore the observations began. Then S'x,2006/'S'j(:,2005 > 20, 
while S'/i?,2006/'S'/fl,2005 > 0.6. The ratio of the submil- 
limeter delays is T2006/2I2005 = 4.8. Assuming that the 
flares are created with similar magnetic field strengths, 
these ratios imply an upper limit on of —2.2. This 
decreases to —2.6 if we assume that the X-ray emission 
is a factor of 2 below the upper limit. Reversing the 
argument, the weakest X-ray flares that can be reliably 
detected correspon d to an excess of app roximately twice 
the quiescent flux. lEckart et al.l ()2004 h observed such a 
flare accompanied by a 4 mJy IR flare (K band). Based 
on the parameters of the 2006 flare and this upper limit 
on we expect that flares with IR-submillimeter de- 
lays smaller than 50 minutes should not show measurable 
X-ray emission. Flares detected in the ongoing coordi- 
nated monitoring campaigns should be able to test this 
relationship in detail. 

A relationship between the ratio of X-ray and IR fluxes 
and the submillimeter delay is expected even if the ex- 
pansion of the flaring region does not follow the power- 
law form assumed above. An example is a plasma re- 
gi on expanding along a je t governed by the equations 
of iFalcke fc Markofil ()2000r ). which account for the accel- 
eration due to the pressure gradient. Because the syn- 
chrotron opacity and the X-ray/IR flux ratio are propor- 
tional to the column density (rigi?) of the plasma, the 
connection is imposed by assumption and the physics of 
the expansion merely determine the form of the correla- 
tion within the limits imposed by magnetic field variabil- 
ity. 

5. SUMMARY 

We have reported the first measurements of a flare of 
Sgr A* observed at submillimeter, IR, and X-ray wave- 
lengths. Interestingly, the submillimeter flare is found 
more than an hour after the X-ray and IR flares. A large 



delay is also found between an IR and submillimeter flare 
(8) in 2005, although the identification of the submillime- 
ter flare with the IR event is less certain. We find the 
spectral and polarization changes in the flare to be con- 
sistent with expansion of a region of ene rgetic plasma, 
althou gh the adiabatic expansion model of lvan der LaanI 
(|1966f ) is not a good fit to this well-observed flare. In- 
dependent of the details of the expansion, this paradigm 
predicts that the delay between the submillimeter and 
short-wavelength flares should be related to the ratio of 
the synchrotron (IR) and SSC (X-ray) luminosities. Such 
a relationship should be testable in the ongoing multi- 
wavelength Sgr A* monitoring campaigns. Campaigns 
including short-wavelength VLBI (1—7 mm) would also 
provide a test of the expansion model, as the long de- 
lays we observe imply expanded so urce sizes compa- 
rable to the measured intrinsic size (|Shen et al.l [20051 : 
iBower et al.ll2006D . The IR and X -ray properties of the 
flares are well modeled by a homogeneous synchrotron- 
SSC source, although incorporation of the submillimeter 
data likely requires a more detailed treatment. Incorpo- 
ration of such expanding flares into the existing static 
models of Sgr A* will be essential for further progress in 
understanding the flares and the accretion region. 



D. P. M. thanks Arieh Konigl, John Raymond, and 
Brant Robertson for enlightening discussions. F. K. B. 
received support for this work from NASA through 
Chandra Award No. G05-6093X and G06-7041X, is- 
sued by the Chandra X-ray Center under contract NAS8- 
03060, and SAO Award No 2834-MIT-SAO-4018. Some 
of these data were obtained at the W.M. Keck Observa- 
tory, which is operated as a scientific partnership among 
the California Institute of Technology, the University of 
California and the National Aeronautics and Space Ad- 
ministration. The Observatory was made possible by 
the generous financial support of the W.M. Keck Foun- 
dation. The National Radio Astronomy Observatory is 
a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, 
Inc. The Submillimeter Array is a joint project be- 
tween the Smithsonian Astrophysical Observatory and 
the Academia Sinica Institute of Astronomy and Astro- 
physics and is funded by the Smithsonian Institution and 
the Academia Sinica. The CSO is supported by the NSF 
under contract AST 05-40882. We thank David Chuss 
and Larry Kirby for assistance with the July 2005 obser- 
vations at CSO. All of the ground-based data presented 
here were obtained from Mauna Kea observatories, a tes- 
tament to the astronomical importance of this site. We 
are grateful to the Hawaiian people for permitting us to 
study the universe from this sacred summit. 

Facilities: CXO (ACIS-I), KeckiH (NIRC2), CSO 
(SHARC-II), SMA (Polarimeter) 



REFERENCES 



Aitken, D. K., Greaves, J., Chrysostomou, A., Jenness, T., Holland, 
W., Hough, J. H., Pierce-Price, D., & Richer, J. 2000, ApJ, 534, 
L173 

Alexander, T. 1997, in Astrophysics and Space Science Library, 
Vol. 218, Astronomical Time Series, ed. D. Maoz, A. Sternberg, 
& E. M. Leibowitz, 163 

Baganoff, F. K. et al. 2001, Nature, 413, 45 



— . 2002, BAAS, 34 
— . 2003, ApJ, 591, 891 

Belanger, G., Goldwurm, A., Melia, F., Ferrando, P., Grosso, N., 
Porquet, D., Warwick, R., & Yusef-Zadeh, F. 2005, ApJ, 635, 
1095 

Belanger, G., Terrier, R., de Jager, O. C., Goldwurm, A., & Melia, 
F. 2006, Journal of Physics Conference Series, 54, 420 



X-ray, IR, and Submillimeter Flare of Sgr A* 



9 



TABLE 1 
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(cm-3) 
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Fig. 1. — Observing windows for the four observatories on 2005 July 31 (left) and 2006 July 17 [right). 
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Fig. 2. — Left: SMA light curve from 2005 July 31. The calibrator was J1733-130 and J1744— 312 has been used as a test source to verify 
the calibration. Right: Light curve from 2006 July 17, with calibrators J1626— 298 and J1924-292, and test source J1733— 130. 
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Fig. 3. — Flux density observed during the 2006 July 17 flare in all three bands. For the CSO 850 fim data the flux density measurements 
are shown at their full temporal resolution {vertical bars) and re-binned into 10 minute averages. The zero point of the 850 fj,m flux density 
scale is uncertain by 1 Jy due to confusion with the surrounding dust emission. The 1.3 mm polarization measured by the SMA is shown 
in the bottom panel. 
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Fig. 4. — Flux density observed on 2005 July 31 in all three bands. The CSO data {bottom) were obtained in three different filters, 
and the frequent switching between wavelengths makes correlation analysis difficult, although the 450 fira data generally follow the trend 
observed at 1.3 mm by the SMA {second from bottom). These have been separated for clarity and plotted to the same scale. The zero point 
of the flux density scales at 850, 450, and 350 fim are uncertain by 1, 0.5, and 1 Jy, respectively, due to confusion. No flare is detected in 
the X-ray observations during this interval, as reported in lHOTl The X-ray and IR data are plotted on the same scale as Figure [3] 
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Fig. 5. — ZDCFs for the available light curves for the 2005 and 2006 flares. Only observations that span the apparent flare peak in each 
band have been cross-correlated. The top panel shows the cross-correlation of the 1.3 mm data (Figure |4] second from bottom) with the 
spectral-average IR light curve. The latter is the combination of the H, K' , and V light curves with the first two scaled by the a = —0.62 
mean spectral index to the L' flux density scale to generate a single light curve with superior sampling. The bottom three panels show the 
three cross-correlations of the 2006 X-ray, 850 /^m, and 1.3 mm data in Figure[3] A positive lag indicates that structure in the second data 
set appears after that in the first. 
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Fig. 6. — Polarization changes during the 2006 July 17 flare at 1.3 mm. The data are binned in half hour intervals (3—6 samples). Top: 
Stokes intensities after subtracting the values in the 630 UT bin (/o=2.98 Jy, Qo=— 89 mjy, C/o=— 128 mjy). The remaining emission is 
ascribed to the flare (the two points that precede the flare are also shown). Bottom: Polarized emission in the flare. Excess polarization is 

and the polarization fraction (m^ 



calculated as Pexcess = \/{Q - Qo)^ + {U 
last points have large and uncertain m, 



excess 



Uq)'^ and the polarization fraction (mexcess) is the ratio of Pj 
; these points and their errors have been scaled by 1/4. 
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Fig. 7. — SED of the peak emission in the 2005 and 2006 flares. The X-ray emission for the 2005 flare is shown as an upper hmit at 
4 keV {triangle). Also shown are three synchrotron-SSC models for the IR and X-ray emission from the flares. The 2006 X-ray flare and 
presumed peak IR emission (chosen to match the 2005 flare maxima) are well fit by model 1 {solid line), while model 2 {long-dashed line) 
approximates the 2006 flare at the start of the IR, coverage. Model 3 { shor t-dashed line) matches the 2005 IR flare and falls below the X-ray 
upper limit. Model parameters are given in Tableland discussed in i|4.4l The submillimeter peaks arc not flt by these models because the 
bulk of the submillimeter photons and low-energy electrons, observed as the quiescent emission, are not accounted for in the flare model, 
and because these peaks occur long after the IR and X-ray peaks. 



